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SUMMARY 


A method  is  presented  for  determining  the  lateral. -stability  deriva- 
tives, transfer-function  coefficients,  and  the  modes  for  lateral  motion 
from  frequency-response  data  for  a rigid  aircraft.  The  method  is  based 
on  the  application  of  the  vector  technique  to  the  equations  of  lateral 
motion,  so  that  the  three  equations  of  lateral  motion  can  be  separated 
into  six  equations.  The  method  of  least  squares  is  then  applied  to  the 
data  for  each  of  these  equations  to  yield  the  coefficients  of  the  equa- 
tions of  lateral  motion  from  which  the  lateral -stability  derivatives  and 
lateral -motion  transfer-function  coefficients  are  computed.  Two  numeri- 
cal examples  are  given  to  demonstrate  the  use  of  the  method. 


INTRODUCTION 


In  the  reduction  and  generalization  of  flight-test  data,  whether  for 
loads,  stability,  or  control  purposes,  the  airplane  stability  derivatives 
and  the  coefficients  of  the  transfer  functions  are  often  required.  A 
great  deal  of  emphasis,  therefore,  has  been  placed  on  the  development  of 
analytical  methods  for  reducing  flight  data  to  obtain  these  basic  deriva- 
tives and  coefficients. 

A number  of  recent  methods,  for  example  references  1 to  k,  are  now 
available  for  analyzing  longitudinal  maneuvers  and  determining  the 
longitudinal-stability  derivatives  and  transfer-function  coefficients 
from  flight  data.  References  1 and  2 present  methods  of  determining  the 
longitudinal-stability  derivatives  and  transfer  functions  directly  from 
transient  data.  Reference  3 reduces  data  for  longitudinal  motion  deter- 
mined from  the  forced-oscillation  technique  by  means  of  circle  diagrams 
to  longitudinal -stability  derivatives  and  frequency  response.  Mueller, 
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in  reference  5,  was  one  of  the  earliest  to  use  vector  representation  in 
the  equations  of  longitudinal  motion  to  represent  the  derivatives  and 
integrals  of  the  variables.  Schumacher,  reference  4,  represented  the 
frequency  responses  to  longitudinal  motion  as  vectors  and  substituted 
them  into  the  equations  of  longitudinal  motion  and  the  transfer  func- 
tions. He  then  applied  the  method  of  least  squares  to  these  vector 
equations,  a method  which  he  found  very  effective  in  determining  certain 
of  the  longitudinal -stability  derivatives  and  transfer-function 
coefficients . 

The  problem  of  analyzing  lateral  motions,  however,  has  not  received 
the  same  amount  of  attention  as  that  for  longitudinal  motion,  perhaps 
because  it  is  more  complicated.  Several  analytical  investigations  have 
been  undertaken  and  a few  methods  have  been  proposed  such  as  the  circle - 
diagram  method  (ref.  6),  the  step-function-response  method  (ref.  6)  in 
which  the  response  of  an  airplane  to  a step  deflection  of  the  rudder  or 
aileron  is  analyzed,  and  the  free -oscillation  method  (ref.  7)  in  which 
the  period  and  damping  of  the  free  vibrations  of  the  aircraft  due  to  a 
pulse -type  input  are  analyzed.  Since  the  usefulness  of  these  methods 
is  limited  by  the  number  of  derivatives  which  can  be  extracted,  there  is 
still  a need  for  a more  general  method  of  analysis  that  will  extract  all 
the  significant  lateral -stability  derivatives  from  flight  data. 

It  is  the  purpose  of  this  paper  to  present  a method  for  determining 
the  lateral -stability  derivatives  of  a rigid  airplane  and  to  illustrate 
its  use  by  applying  it  to  two  examples.  The  method  is  based  on  the 
vector  representation  of  the  frequency  responses  to  lateral  motions. 

This  vector  approach  permits  separation  of  each  of  the  equations  of 
lateral  motion  into  a real  and  imaginary  equation.  A least-squares 
method  is  then  applied  to  the  data  in  each  of  these  equations  or  combi- 
nations thereof  to  yield  the  coefficients  of  the  equations  of  lateral 
motion.  The  lateral-stability  derivatives  and  transfer-function  coef- 
ficients are  then  determined  from  these  coefficients  and  the  known  air- 
craft mass  parameters . 

The  method  is  applied  to  two  specific  examples,  one  in  which  the 
frequency  responses  to  a rudder  input  are  known  and  one  in  which  the 
transient  responses  to  aileron  deflection  are  known.  In  the  latter  case, 
the  frequency  responses  were  obtained  from  the  transient  motions  by  two 
methods  and  the  stability  derivatives  computed. 

An  attempt  has  been  made  to  schedule  the  procedure  so  as  to  reduce 
the  dependence  of  the  results  obtained  from  this  method  on  the  deriva- 
tives that  can  be  least  accurately  obtained  from  the  particular  data 
being  analyzed;  however,  further  improvements  may  be  made  as  further 
experience  is  gained  in  the  application  of  the  method. 
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®y 

Ap,  A^,  A^,  Ag^ 
BP>  B<p>  Bay 

b 


CL 


Cn 

Cz 


Cl,  c2,  c^, 


D 


Fl>  F2>  f3 
Gg,  Gj 


g 


I 


X 


IZ 


SYMBOLS 

lateral  acceleration,  ft/sec2 
parameters  defined  by  equations  (8) 

parameters  defined  by  equations  (9) 

wing  span,  ft 

lift  coefficient,  L/qS 

yawing -moment  coefficient,  N/qSb 

rolling -moment  coefficient,  L'/qSb 

lateral -force  coefficient,  Lateral  force/qS 

coefficient  of  transfer  functions 
(defined  in  table  II) 

differential  operator,  -£7- 

dt 

natural  logarithmic  base 

forcing-function  coefficients  representing  rudder 
effectiveness  (defined  in  table  I) 

forcing-function  coefficients  representing  aileron 
effectiveness  (defined  in  table  I) 

acceleration  due  to  gravity,  ft/ sec2 

moment  of  inertia  about  stability  X-axis,  mk^2, 
slug-ft2 

moment  of  inertia  about  stability  Z-axis,  mkz2, 
slug-ft2 


} 
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Xv-  product  of  inertia  referred  to  stability  axes  (nega- 

tive when  the  positive  direction  of  the  X principal 
axis  is  inclined  above  the  flight  path,  i.e.,  when 
t]  is  positive) 

i = 4-1 

K^,  K2,  Kj,  . . . stability-derivative  coefficients  of  the  equations 

of  lateral  motion  (defined  in  table  I) 


KZ 


Kxz 


krr 

^0 

L 

L' 


M 


m 

N 

p or  0 


radius  of  gyration  about  stability  X-axis,  ft 
radius  of  gyration  about  stability  Z-axis,  ft 
nondimens ional  radius  of  gyration  about  longitudinal 


stability  axis,  (ky^jbj  ® cos^T)  + ^kz  I bj  ^ sin^rj 


nondimens ional  radius  of  gyration  about  vertical 
stability  axis, 


(kz0/b)  2 cos2t>  + (kxG/b)  2 sin2li 


nondimens ional  product-of -inertia  parameter, 


- (%/»)  2 - (Kx0/1 ) 2 


cos  t]  sin  T) 


radius  of  gyration  about  principal  longitudinal 
axis,  ft 

radius  of  gyration  about  principal  vertical 
axis,  ft 


lift,  lb 

rolling  moment,  ft -lb 
Mach  number 

aircraft  mass,  W/g,  slugs 
yawing  moment,  ft -lb 

incremental  rolling  angular  velocity  about 
X-axis,  radians/sec 
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CL 

r or  if 

R 

S 

s 

t 

v 


V 

W 


X,  Y,  Z 


a 

3 


6r 


^b 

P 

T 

$ 

* 

a) 


dynamic  pressure,  |pV2,  lb/sq  ft 

incremental  yawing  angular  velocity  about 
Z-axis , radians/sec 

amplitude  ratio 

wing  area,  sq  ft 

Laplace  transform  variable 

time,  sec 

incremental  component  of  velocity  V alone 
the  Y-axis,  ft/sec 

true  airspeed,  ft/sec 

aircraft  weight,  lb 

airplane  stability  axes  (see  fig.  1) 

angle  of  attack,  radians 

angle  of  sideslip,  v/v,  radians 

aileron  control  deflection,  radians 

rudder  control  deflection,  radians 

inclination  of  principal  longitudinal  axis  of  inertia 
with  respect  to  flight  path  (positive  when  the  posi- 
tive direction  of  the  X principal  axis  is  inclined 
above  the  flight  path) 

nondimens ional  mass  parameter  used  for  lateral 
equations,  m/pSb 

air  density,  slugs/cu  ft 

time  parameter,  m/pSV,  sec 

phase  angle , radians 

angle  of  roll,  radians 

angle  of  yaw,  radians 

angular  frequency,  radians/sec 
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Bar  notation: 


a bar  over  letter  represents  maximum  value 

| a | bars  on  sides  of  symbol  represent  absolute  value 

Matrix  notation: 


[] 

{} 

The 
tion  as, 


square  matrix 
column  matrix 

lateral-stability  derivatives  are  expressed  by  subscript  nota- 

dC7  dCL,  dCn 

for  example:  C7  = C.,  = — — , and  Cn  = — -. 

z3  ^3  "p  ^b  ^ ^rb 

2V  °2V 


Phase  angles  are  also  indicated  by  subscript  notation  as  $0  is 

p6r 

the  phase  angle  between  the  input  rudder  deflection  and  the  output  side- 
slip angle. 


DEVELOPMENT  OF  METHOD 
Equations  of  Lateral  Motion 


The  equations  of  lateral  motion  based  on  the  stability  axes  (as  shown 
in  fig.  l)  and  on  the  usual  assumptions  of  linearity,  small  angles,  and 
maneuvers  which  start  from  a level-flight  condition  are  usually  written 
in  the  following  form: 

(a)  For  the  sideslipping  motion, 

(D  + Kx)3  - K20  + D+  = Fjb^t)  (1) 

or,  since  = v(D3  + Dt  - K20), 

iay  + K!3  = F!&r(t)  (la) 

(b)  For  the  rolling  motion, 

K53  + (d2  + K4d)  0 - (K5D2  + K6d)  t = F25r(t)  + G26a(t)  (2) 
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(c)  For  the  yawing  motion, 

-y  - (KqD2  + K^d)  ( + (d2  + K10Dj  * = F55r(t)  + G^Jt)  (3) 


where  the  K,  F,  and  G coefficients  are  defined  in  table  I.  The 
forcing  functions  of  these  equations  are  written  for  rudder  inputs  and 
aileron  inputs  combined;  for  the  sake  of  simplicity  in  developing  and 
demonstrating  the  method,  only  rudder  inputs  F^6r(t),  Fg5r(t), 

and  F^5r(t)  are  considered.  The  method,  however,  applies  equally  well 

to  both  aileron  and  rudder  inputs  as  is  shown  subsequently. 


Vector  Interpretation  of  Frequency  Response 

In  many  current  studies  of  airplane  dynamics,  the  frequency  response 
of  the  aircraft  is  determined  from  flight  data.  Having  the  data  avail- 
able in  the  frequency  plane  offers  certain  advantages  over  having  the 
data  in  the  time  plane,  since  use  can  then  be  made  of  the  vector  inter- 
pretation of  frequency  response.  The  development  which  follows  will 
assume  that  the  data  to  be  analyzed  are  available  in  frequency-response 
form.  In  the  event  that  only  transient-response  data  are  available,  the 
transformation  to  the  frequency  plane  can  be  made  from  a selection  of  one 
of  the  several  methods  compared  in  reference  8. 

The  vector  technique  is  applied  to  the  determination  of  lateral- 
stability  derivatives  in  the  following  manner.  If  an  input  5r  to  a 
linear  system  (a  system  described  by  a linear  differential  equation)  is 
considered  to  have  a sinusoidal  variation  of  frequency  u>,  it  can  be  con- 
sidered to  be  a vector  of  unit  magnitude  lying  along  the  real  axis  of  the 
complex  plane  as  shown  in  figure  2(a).  The  first  derivative  of  the  sinus- 
oidal input  D5r  is  then  obtained  by  multiplying  the  amplitude  by  w and 
rotating  the  resulting  vector  90°  counterclockwise  in  the  complex  plane 
which  is  equivalent  to  advancing  the  phase  angle  by  9O0.  Each  succeeding 
higher  derivative  is  found  by  multiplying  by  o>  and  rotating  the  vector 
90°  in  the  complex  plane;  thus, 


= 1 

= iu> 

(5) 

= -o£ 

(6) 

• i 

( 


I 

\ 


K 


\ 
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The  steady-state  output  or  response  of  the  linear  system  to  a sinus- 
oidal input  may  be  considered  to  be  a vector  of  magnitude  equal  to  an 
amplitude  ratio  R and  having  a direction  angularly  displaced  from  the 
real  axis  by  a phase  angle  0.  Thus,  in  the  complex  plane,  the  vector 
for  sideslip  may  be  represented  as 


3 = 


i_ 

‘V 

e 1 = 

3 

&r 

5r 

cos 


°05r  + 1 sin<^5r)  = + iBP 


(7) 


By  definition, 

Ap  = 

A^  = 


cos  Of 


&r 


cos  0 


B3" 


5r 

J- 

6r 


sin 


sin  0^ 


sL 


cos  0^ 


> (8) 


A-  s 


B0" 


t 

6r 


cos  0, 


Ay&T 


sin  0 


B»  = 

®y 


h. 

&r 


sin  0a. 


y5r 


The  derivatives  of  3 are  represented  as 


(9) 


Dp  = 0) 


= CL) 


ifXr+2 


-sin  0^^  + i cos  0p^ 


= -aiBp  + iuAp 


(10) 
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and. 


Similarly,  for  roll  (see  fig.  2(b)),  the  equations  corresponding  to  equa- 
tions (7),  (10),  and  (ll)  would  be 


for  yaw, 


and,  for 


0 = + iB0  (12) 

wf)  = -oiB0  + ioA^  (13) 

D 20  = -o>2A^  - io>2B0  (14) 

\|f  = A^  + iBy  (15) 

IN'  = -o%  + iojfy  (16) 

D2^  = - io^  (17) 

lateral,  acceleration, 

®y  = Aay  + iBay  (18) 


Substituting  equations  (4)  to  (18)  into  equations  (l),  (la),  (2),  and  (3) 
and  equating  the  real  and  imaginary  values  in  each  equation  yields  the 
following  eight  simplified  equations : 


i 


\ 


K 
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From  equation  (l): 


KjBp  - KgB^  = -w(Ap  + Aj)  (19) 

and 

KxAp  - = Fi  + o>(Bp  + B^)  (20) 

From  equation  (la): 

VKjBp  = -B^  (21) 

VK^  = VFL  - A^  (22) 

From  equation  (2): 

K^Bp  + K4aVy  + ICpU)2^  - K^cuA^  = a^B^  (23) 

and 

K3A0  - K4o£0  + K5^  + K6mBt  - F2  = o)2^  (24) 

From  equation  (3): 

-KyBp  + Kgo^BjjJ  - KytiA^  + K^quA^  = (23) 

and 

-KyAp  + Kqo)2A^  + KgoB^  - K^qoB^  - Fj  = cd2A^  (26) 


In  equations  (19)  to  (26)  the  Ap^^^  and  Bp,0,\jr,ay  terms  wil1 

be  available  at  the  particular  values  of 
curves.  It  will  be  assumed  that  the  K2 
velocity  measurement  and  the  K5  and  K8 


a)  from  the  frequency-response 
term  is  determined  from  the 
terms,  which  are  equal  to  the 


ratios  Iyz/^x  and  IXz/IZ'  respectively,  are  known  from  either  weight 


and  balance  calculations  or  measurements  so  that  only  the  values 
of  ^3^6,7,9,10  and  ^1,2,3  are  to  ^ determined  from  the  data. 

(it  should  be  noted  that  even  though  and  Kq  are  assumed  to  be  known 


a need  still  exists  for  a simple  method  of  determining  aircraft  moments 
of  inertia  and  the  location  of  the  principal  longitudinal  axis.) 


{ 
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Equations  (19)  to  (2 6)  are  now  used  to  compute  the  K and  F coef- 
ficients. The  advantage  of  separating  each  of  the  equations  of  lateral 
motion  into  real  and  imaginary  equations  now  becomes  apparent.  The  two 

derivatives  most  difficult  to  determine  are  C7  and  Cn  contained  in 

lr  “p 

the  Kg  and  Kg  coefficients,  respectively.  The  Kg  coefficient  may- 


be eliminated  between  equations  (23)  and  (24);  similarly  the  Kg  coef- 
ficient may  be  eliminated  between  equations  (25)  and  (26).  The  equations 
thus  obtained,  when  used  to  compute  the  unknown  coefficients,  result  in 
better  conditioned  matrices  and  more  accurate  results.  In  the  usual  case 
it  has  been  found  unnecessary  to  perform  the  elimination  of  the 
K^  coefficient. 


As  a first  step  in  solving  for  the  K and  F coefficients,  equa- 
tions (19)  to  (26)  are  fitted  to  the  frequency-response  data  over  a range 
of  values  of  a)  by  the  method  of  least  squares.  The  theory  of  least 
squares  is  well-known  and  is  derived  in  many  textbooks  (e.g.,  ref.  9)  and 
no  attempt  is  made  here  to  repeat  its  derivation. 


The  application  of  the  least-squares  method  to  equations  (19)  to  (26) 
converts  these  equations  to  their  computational  form.  It  is  recommended 
that,  if  the  lateral -acceleration  frequency  response  is  available,  equa- 
tions (21)  and  (22)  be  used  in  preference  to  equations  (19)  and  (20)  for 
the  computation  of  Kp  and  F^. 

Equation  (19)  may  be  expressed  simply  as 


and  yields  Kj_.  Similarly,  equation  (20)  can  be  expressed  as 


^(KiAp  - KgA^j  - 


o£p 


(27) 


(28) 


and  yields  F-^. 
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For  the  same  variables  K]_  and  Fj_,  the  least-squares  forms  of 
equations  (21)  and  (22)  for  use  with  the  lateral -acceleration  frequency 
response  are: 


* ’*£.&*%  <50) 

Equation  (25)  can  be  expressed  as 


V 

^ '-N 

p (A  - 

A 

... 

^ (A -■%“%) 

2 f-Bp)  ,'«**)  , 
J3.  J J 

K10 

(A  - ^ A) , K) , 

' j jj 

and  yields  Kj,  K^,  and  K^g. 

Equation  (26)  can  be  expressed  as 

r5  - i J_  K*s  + + lt9“®d  - - “%) , (52) 

J=o>L  v « 

and  yields  Fj . 

In  order  to  evaluate  the  remaining  constants,  K^,  K^,  Kg,  and  Fg, 

the  coefficient  Kg  is  first  eliminated  between  equations  (23)  and  (24). 
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This  step  is  taken  in  order  to  obtain  a more  accurate  value  of  the  coef- 
ficient K5,  which  is  less  significant  tinder  these  conditions,  and  to 
provide  a better  conditioned  matrix  of  coefficients  in  the  computation 
of  Kj,  K4,  and  Fg«  Upon  elimination  of  the  Kg  term  between  equa- 
tions (23)  and  (2k),  the  following  relation  is  obtained: 

- 2 

+ ApA^  + K^oj^A^B^  - B^A^  - FgA^  = + A^A^  - KtjU^ 


Upon  the  application  of  the  method  of  least  squares,  equation  (33) 
becomes 


^ j^e8*  ♦ V*)  ] 2 - “***♦)  j + V*) } ^ ('**)  j (V*  + V*)  3 Ki 

p.  M ♦ V*)  j “¥*♦  - «***♦)  j ^ - “V*/1]2  j(“V*  ' “V*;  j,Kl* 


- S‘B8|^  j(V*  * V*)  j 
♦ <Ay\  - ■ “*A)j 


+ “ Y»  - v H LWj 


which  yields  K^,  K^,  and  F2. 

The  Kg  coefficient  may  now  be  obtained  from  equation  (24)  which, 
after  the  method  of  least  squares  has  been  applied,  may  be  expressed  as : 

y*  (f2  + cii2Ad  - K^aj2A+,  + %a£d  - KjApl  (uB|) 

1%  = ±3d - <3?: 

kM1 
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An  alternate  expression  for  Kg 
which  may  he  expressed  as 


may  he  obtained  from  equation  (23) 


Values  of  Kg,  which  define  the  derivative  C j , are  important  in 

the  spiral  mode.  In  order  to  obtain  accurate  values  of  Kg,  aircraft 

motions  which  bring  out  the  full  effect  of  the  long-period  spiral  mode 
should  be  analyzed  to  obtain  accurate  frequency-response  data  near  cs  = 0. 

In  the  case  where  the  solution  for  C„  is  found  to  involve  in- 
i’ 

conditioned  matrices,  should  be  eliminated  between  the  real  and 

imaginary  equations  (25)  and  (26). 

With  all  the  K and  F coefficients  of  the  equations  of  motion 
determined,  the  lateral-stability  derivatives  can  be  computed  from  the 
K and  F coefficients  and  the  known  aerodynamic  parameters  by  use  of 
the  definitions  of  table  I.  The  transfer-function  coefficients  can  be 
calculated  from  the  definitions  in  table  II  and  the  transfer  functions 
and  modes  of  lateral  motion  can  be  obtained  from  the  equations  in 
appendix  A. 


SUGGESTED  PROCEDURE 


In  order  to  aid  in  the  application  of  this  method  to  the  analysis 
of  flight  data,  a suggested  step-by-step  procedure  which  is  presented  in 
this  section  has  been  worked  out.  An  effort  has  been  made  to  schedule 
the  procedure  so  as  to  reduce  the  dependence  of  the  results  on  the  deriva- 
tives that  can  be  least  accurately  obtained  from  the  particular  data  being 
analyzed.  Alternate  steps  are  suggested  where  it  was  found  that  some  par- 
ticular derivative  might  be  more  accurately  determined  by  one  or  the  other 
of  two  approaches  under  certain  conditions.  No  weighting  of  specific 
groups  of  data  is  employed  in  the  least-squares  procedure;  weighting  can, 
however,  be  employed  when  it  is  considered  desirable  to  put  more  dependence 
on  data  regarded  as  more  reliable. 
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As  a further  demonstration  of  the  method,  two  numerical  examples 
have  been  carried  out  according  to  this  procedure.  In  one  example,  pre- 
sented in  appendix  B,  the  lateral-stability  derivatives  are  calculated  for 
an  airplane  whose  transient  responses  to  aileron  deflection  are  assumed 
to  be  known.  This  example  was  chosen  in  order  to  gain  some  insight  into 
the  dependence  of  the  method  on  the  accuracy  with  which  the  frequency 
response  can  be  obtained  from  the  transient  response. 


The  other  example  is  for  the  hypothetical  rigid  airplane  whose  mass 
and  geometric  parameters  are  listed  in  table  III.  The  assumed  frequency- 
response  data  of  sideslip  angle  3,  roll  angle  0,  yaw  angle  +,  and 
lateral  acceleration  a y for  a rudder-deflection  input  such  as  might  be 
obtained  from  analysis  of  flight  data  are  listed  in  table  IV  and  plotted 
in  figures  3 to  6.  The  analysis  of  these  data  has  been  used  to  illustrate 
the  following  suggested  procedure. 


(l)  Tabulate  parameters  and  working  equations.-  Tabulate  the  air- 
plane parameters  and  any  stability  derivatives  that  are  known  or  can  be 
computed  directly  from  those  that  are  known  as  illustrated  in  table  III. 
Tabulate  the  least-squares  equations  (27)  to  (32)  and  (34)  to  (36)  with 
any  known  terms  on  the  right-hand  side  as  has  been  done  in  equations  (2'Ja) 
to  (32a)  and  (34a)  to  (36a)  in  table  V for  the  example.  (Since  CUp  was 

taken  to  be  zero  in  the  example,  all  terms  containing  were  dropped; 

however,  in  the  example  of  appendix  B,  was  included  and  the 

K9  term  was  determined.) 


(2)  Tabulate  frequency-response  data.-  Tabulate  the  amplitudes  and 
phase  angles  of  the  four  lateral  variables  at  the  values  of  o>  to  be 
used  in  the  analysis.  For  simplicity  in  the  example,  as  shown  in  table  IV, 
10  integral  values  of  cd,  evenly  distributed  over  the  range  of  the  data, 
were  chosen;  however,  in  cases  where  more  accurate  points  occur  at  non-, 
integral  values  of  ou,  it  Is  advisable  to  use  such  points,  where  possible. 
The  range  of  values  of  cd  should  be  restricted  to  the  rigid  response, 


but  the  point  cu 


0 


should  be  avoided  because 


is  infinite. 


(3)  Compute  vector  components  of  data.-  Compute  the  vector  compo- 
nents from  the  data  listed  in  step  2 and  tabulate  them  as  shown  in 
table  IV.  It  is  advisable  to  use  at  least  5-place  tables  of  sine  and 
cosine  values  in  this  computation.  Additional  columns  of  oAp,  cdBp, 

<uA^,  . . . and  o^A^,  o^B p,  u£a^,  . . . will  be  useful  in  subsequent 
operations  but  are  not  listed  for  the  example. 


\ 
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(4)  Compute  Kp.-  Compute  Kp  from  equation  (27a).  (See 

table  Vl(a).)  In  this  step  and  in  subsequent  steps,  summation  of  the 
corresponding  elements  of  the  indicated  products  of  the  data  columns 
evaluates  the  matrices  of  coefficients  and  knowns . The  resulting  simul- 
taneous equations  are  readily  solved  by  any  standard  method. 

(4a)  Alternate  step  4.-  Compute  Kp  from  equation  (29a).  (See 
table  VI ( c ) . ) 

(3)  Compute  F-[.~  Compute  Fp  from  equation  (28a)  with  the  values 
of  Kq  from  step  (4).  (See  table  Vl(b).) 

(3a)  Alternate  step  3.-  Compute  Fp  from  equation  (30a)  with  the 
values  of  Kp  from  step  (4a).  (See  table  Vl(d).) 

(6)  Compute  Kfj  and  KpQ.-  Compute  K ^ and  Kpg  from  equation  (31a). 
(See  table  Vl(e).) 

(7)  Compute  F^. - Compute  from  equation  (32a)  with  the  values 

of  K 7 and  Kp0  from  step  (6).  (See  table  Vl(f).) 

(8)  Compute  K^,  Kq,  and  F2.-  Compute  K^,  Kq.,  and  F2  from 

equation  (34a).  Because  of  the  combination  of  two  equations  to  obtain 
equation  (34)  in  the  method,  the  combination  of  data  columns  is  more 
complicated  in  this  step  than  in  others.  (See  table  VI (g).) 

(9)  Compute  Kg.  - Compute  Kg  by  one  of  the  following  methods: 

(a)  If  Fg  has  not  been  eliminated,  calculate  Kg  from  equa- 
tion (35a)  with  the  values  of  Kj,  K^,  and  F2  from  step  (8).  (See 
table  Vl(h) . ) 

(b)  If  F2  has  been  eliminated  calculate  Kg  from  equation  (36a). 

(See  table  Vl(i).) 

(10)  Compute  stability  derivatives.-  Compute  the  lateral-stability 
derivatives  from  the  relationships  in  table  I and  the  now  complete  set  of 
values  of  K and  F.  (The  computed  and  known  values  of  the  K 

and  F coefficients  and  the  stability  derivatives  for  the  example  are 
given  in  tables  VII  and  VIII,  respectively.) 


I 
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(11)  Compute  transfer-function  coefficients.-  Compute  the  transfer- 
function  coefficients  from  the  relationships  given  in  table  II  and  the 
values  of  K and  F.  (The  computed  and  known  values  of  the  transfer- 
function  coefficients  for  the  example  are  given  in  table  IX.) 

(12)  Compute  the  modes.-  The  modes  of  lateral  motion  may  now  be 
computed,  if  desired,  by  use  of  the  equations  in  appendix  A.  (The  values 
of  the  lateral  modes  for  the  example  are  given  in  table  X.) 

(13)  Accuracy  check.-  When  the  time  derivatives  are  unknown,  the 

check  on  the  computations  is  obtained  in  the  following  manner:  Calculate 

the  frequency  response  from  the  equations  in  appendix  A and  the  computed 
coefficients  of  step  (ll).  Compare  these  values  with  the  data  of 

step  (2). 


DISCUSSION 


Data  requirements.-  The  method  of  this  paper  is  based  on  the  assump- 
tion that  equations  (1),  (2),  and  (3)  represent  the  lateral  motions  of  the 
airplane . These  linear  differential  equations  are  satisfied  by  small 
variations  in  roll,  yaw,  and  sideslip  but  will  not  necessarily  hold  for 
large  changes  in  these  variables  or  in  regions  of  flight  where  nonlineari- 
ties occur.  It  is  therefore  desirable  to  obtain  data  for  transformation 
to  the  frequency  plane  from  small  motions  recorded  within  ranges  of  Mach 
number  and  other  variables  in  which  the  coefficients  are  constant;  also, 
the  range  of  frequency-response  data  to  be  used  must  be  limited  to  values 
for  which  the  aircraft  acts  as  a rigid  body  and  higher  frequency  instru- 
ment inaccuracies  are  avoided.  Thus,  the  assumptions  of  the  method  will 
be  most  nearly  in  accord  with  the  true  physical  nature  of  the  problem. 

In  any  case,  experience  has  shown  that  the  derivatives  obtained  by  this 
method  of  analysis  will  accurately  reproduce  the  frequency  response  frem 
which  they  were  extracted. 

Since  the  present  method  requires  flight  data  in  frequency-response 
form,  either  the  forced-oscillation  technique  or  the  technique  of  oscil- 
lating the  rudder  or  aileron  at  various  frequencies  and  measuring  the 
steady-state  response  in  3,  D 0,  Ity,  and  Sy  must  be  employed  or  the 

transformation  of  a transient  response  from  the  time  plane  to  the  fre- 
quency plane  must  be  made  by  an  appropriate  method.  A comprehensive  sur- 
vey of  methods  for  effecting  this  transformation  from  time  to  frequency 
plane  is  given  in  reference  8 which  evaluates  the  Fourier  integral  method 
as  well  as  the  use  of  devices  such  as  the  Fourier  synthesizer,  Coradi 
harmonic  analyzer,  and  IBM  machines  for  obtaining  frequency  response  from 
transient-response  data. 
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If,  as  is  usually  the  case,  the  transient  responses  in  roll  and  yaw 
are  measured  in  the  form  of  time  histories  of  D0  and  D>|r  rather  than 
of  ^ and  \|r,  it  is  recommended  that  these  time  histories  of  D <f> 
and  Dty  be  converted  to  frequency  response  initially  and  then  the  fre- 
quency responses  of  (f>  and  \|r  be  computed  from  the  frequency  responses 
of  and  Ity  by  use  of  the  relations 


= 0) 

I 

= CL) 

± 

Sr 

5r 

5r 

5r 

X * % * 90°  x ■ V + 900 


Flight  data  are  usually  measured  about  the  body  axes;  before  they 
can  be  used,  however,  either  they  must  be  converted  to  a form  corre- 
sponding to  the  stability  axes  or  the  development  must  be  restated  in 
terms  of  body  axes.  The  following  equations  are  presented  to  be  used 
in  transforming  to  the  stability  axes: 

<f>  = cos  a + sin  a 

• • \ 

i|r  = cos  a - 55^  sin  a 


where  a is  the  angle  of  attack  (angle  between  the  body  X-axis  and  the 
stability  X-axis)  and  the  subscript  b refers  to  the  body  axes.  In 
addition,  the  sideslip-vane  measurements  can  be  converted  to  the  corre- 
sponding values  at  the  center  of  gravity  by  means  of  the  formula 


3 


cos  a 


where  the  subscript  v refers  to  the  vane  and  Z^  is  the  distance  from 
the  center  of  gravity  to  the  sideslip  vane  measured  along  the  body  X-axis. 

If  the  lateral  acceleration  ay  as  recorded  by  an  accelerometer  is 

used,  it  must  be  corrected  for  any  components  of  angular  velocity  and 
acceleration  present  in  such  records  because  of  the  displacement  of  the 
accelerometer  from  the  airplane  center  of  gravity. 
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Sensitivity  of  the  method  and,  possible  simplifications.-  Even  though 
the  foregoing  general  rules  are  observed,  the  question  arises  as  to  how 
the  derivatives  obtained  compare  with  the  actual  derivatives  of  the  air- 
plane, and  how  sensitive  these  derivatives  are  to  inaccuracies  in  the 
frequency-response  data.  These  inaccuracies  stem  from  several  sources 
such  as  measurement  of  the  transient  motions,  transformation  from  time  to 
frequency  plane,  and  reading  errors.  Care  should  be  exercised  to  keep 
these  inaccuracies  to  a minimum. 


No  attempt  is  made  to  present  a comprehensive  error  analysis  in  this 
paper,  but  several  general  remarks  appear  to  be  in  order.  A comparison 
of  the  values  given  in  tables  VII,  VIII,  and  IX  for  the  coefficients  of 
the  equations  of  lateral  motion,  the  lateral-stability  derivatives,  and 
the  transfer-function  coefficients  indicates  essentially  no  errors  due 
to  the  method  itself.  It  has  been  found,  however,  that  the  derived  values 


of  Cv  and  C 
YP 


YSr 


ties  An 


Bn 


are  sensitive  to  slight  variations  in  the  quanti- 

(see  eqs.  (27)  and  (28)).  This  sensitivity 


and  B>t 


^r 

CU1U- 

is  due  to  the  fact  that  the  airplane  angle  of  yaw  is  very  nearly  equal 
to  the  negative  of  the  sideslip  angle  and  hence  accurate  determination 
of  the  lateral  acceleration  from  measurements  of  ijr  and  p is  diffi- 
cult. This  difficulty  can  be  obviated  by  a direct  measurement  of  the 
lateral  acceleration  with  an  accelerometer  and  more  accurate  values 
determined  for  Cy  and  Cy  by  means  of  equations  (29)  and  (30). 


As  brought  out  in  appendix  B,  the  accuracy  with  which  the  values  of 
the  derivatives  Cjr,  Cnr,  and  Cy^  can  be  determined  from  frequency- 

response  data  appears  to  be  dependent  upon  the  inclusion  of  data  for  fre- 
quencies near  the  natural  frequency  of  the  airplane.  This  fact  is  brought 
out  further  by  seme  unpublished  results  that  indicate  that  derivatives 
such  as  Cnr  and  Cy^  are  very  sensitive  to  variations  of  the  frequency 

responses  near  the  airplane  natural  frequency.  In  the  case  of  Cnr  the 


conclusion  is  obvious  since  this  derivative  is  known  to  contribute  largely 
to  the  airplane  Dutch  roll  damping  and  it  is  also  known  that  the  peak  value 


of  the  frequency-response  amplitude  ratio,  say 


is  highly  dependent 


upon  the  airplane  Dutch  roll  damping.  It  appears  therefore  that  an 
accurate  representation  of  the  frequency  responses  at  frequencies  near 
the  airplane  natural  frequency  (for  instance  ±50  percent)  is  probably 
very  important  in  the  determination  of  the  lateral -stability  derivatives. 


In  general,  when  the  number  of  unknown  derivatives  can  be  reduced, 
the  remaining  derivatives  can  be  determined  with  greater  reliability. 

For  this  reason,  several  possible  simplifications  are  presented  for  the 
calculation  of  the  derivatives  C_  , C„  , C,  , and  C,  . The  airplane 

V n8r  V l8a 
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natural  frequency  is,  for  small  values  of  the  product  of  inertia,  well 
approximated  by  the  expression 


p 

off  «>  C 


qSb 


^Z 


If  the  natural  frequency  ajj  is  apparent  from  the  available  transients 
and  the  inertia  parameter  is  known,  a very  good  value  of  Cn  can  gen- 


erally be  obtained  from  this  expression. 


Extrapolation  of  the  amplitude  ratio 


± 


■0 


to  zero  frequency  may 


afford  a means  of  obtaining  the  derivative  Cn&  subsequent  to  deter- 
mination of  Cnp  from  the  preceding  frequency  relation.  The  value  of 

■JU  at  a)  = 0 obtained  by  this  extrapolation  is  not  indicative  of  the 
I °r  I 

actual  static  sensitivity,  but  rather  of  the  apparent  static  sensitivity 

parent  stei 

from  which 


which  includes  the  very  lightly  damped  spiral  mode 
state  response  is  approximately  equal  to  the  ratio 

is  known. 


Cn&  can  be  estimated  if  Cnp 


This  apparent  steady- 

Cn&r 

S ’ 


The  transfer  function  which  relates  the  rolling  velocity  to  aileron 
deflection  is  well  approximated  (particularly  at  low  frequencies)  by  the 
expression 


L 


-i  tan”1  ^ 


(io>) 


Ki.  + io) 


Therefore,  (co)  = -tan"1  — from  which  K 4,  and  subsequently  C 1 , 

Ki.  P 


k4 


(K42  + (lf~J 


1/2 


k4 

can  be  determined.  Also,  by  extrapolating  the  amplitude  ratio 


1 


to 


a)  = 0,  the  parameter  G2  and  hence  C,  can  be  obtained  once  has 

been  determined.  As  was  the  case  for  )■£- 


yield  the  true  value  of 


li 


, this  extrapolation  does  not 


at  u)  = 0,  which  is  actually  zero,  but  gives 


the  apparent  value  due  to  the  lightly  damped  spiral  mode. 
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CONCLUDING  REMARKS 


A method  has  been  presented  for  extracting  lateral-stability  deriva- 
tives from  frequency-response  data  which  have  been  derived  from  aircraft 
transient  responses  to  arbitrary  control  inputs.  In  order  to  demonstrate 
the  use  of  the  method,  the  lateral-stability  derivatives  have  been  calcu- 
lated for  two  hypothetical  airplanes  for  which  the  frequency-response  data 
or  transient-response  flight  data  were  assumed  to  be  known.  Simplifica- 
tions were  proposed  for  obtaining  certain  of  the  derivatives,  and  although 
no  error  analysis  is  presented  in  this  paper,  seme  general  observations 
can  be  made  concerning  the  sensitivity  of  the  method  to  inaccuracies  in 
the  original  data. 

In  view  of  the  limited  experience  in  the  determination  of  lateral - 
stability  derivatives  from  flight  data,  and  particularly  in  using  the 
present  method,  it  appears  that  an  investigation  should  be  made  to  deter- 
mine the  effects  of  some  of  the  errors  inherent  in  analyses  of  transient 
responses,  such  as  recording  accuracy,  reading  accuracy,  and  so  forth, 
on  the  derived  frequency  responses  and  hence  on  the  computed  lateral- 
stability  derivatives.  It  is  noted  here  that,  although  the  methods  pre- 
sented in  this  paper  for  the  extraction  of  lateral-stability  derivatives 
from  frequency-response  data  are  theoretically  correct  and  mathematically 
sound,  the  accuracy  with  which  the  lateral -stability  derivatives  of  an 
airplane  can  be  determined  by  this  method  is  directly  dependent  upon  the 
accuracy  with  which  the  frequency  responses  of  the  airplane  are  known. 
Also,  seme  consideration  should  be  given  to  the  types  of  inputs  likely 
to  afford  good  frequency-response  data  and  to  testing  procedures  whereby 
some  of  the  stability  derivatives  that  axe  most  difficult  to  determine 
may  be  accurately  obtained.  In  addition,  accurate  frequency-response 
data  for  specific  ranges  of  frequency  may  yield  more  accurately  certain 
of  the  stability  derivatives  and  this  possibility  should  be  investigated. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  February  16,  195^* 
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APPENDIX  A 

DEVELOPMENT  OF  EQUATIONS  FOR  TRANSFER  FUNCTIONS,  MODES, 
AND  FREQUENCY  RESPONSE  FOR  LATERAL  MOTION 


In  this  appendix  the  transfer  functions  are  derived  from  the  equa- 
tions of  lateral  motion  and  the  modes  and  frequency  response  axe  expressed 
in  terms  of  the  transfer-function  coefficients. 


Equations  of  Motion 

The  motions  of  a rigid  aircraft  resulting  from  a rudder-deflection 
input  are  assumed  to  be  expressed  by  the  following  three  standard  line- 
arized equations  of  lateral  motion: 

(D  + K^p  - + D|r  = F15r(t) 

(Al) 

+ (1^  + K^D)0  -(i^D2  + KgD)*  = F2&r(t) 

(A2) 

-KyP 

-(kQD2  + +(l£  + K1qd)*  = F5Br(t) 

(A3) 

where  the  K and  F coefficients  are  defined  in  table  I.  The  axes  and 
the  sign  conventions  employed  are  shown  in  figure  1.  On  applying  the 
Laplace  transformation 


F(s) 


f(t)e’stdt 


(A4) 


(where  s = a + io>)  to  both  sides  of  equations  (Al),  (A2),  and  (A3)  and 
assuming  the  initial  conditions  to  be  zero,  there  is  obtained 
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8 + Ki  -K2 

K3  s2  + K^s 

-Ky  -Kgs2  - K^s 


s 

3(s) 

VlM®) 

-y2  - y 

0(s) 

F25r(s) > 

s2  + K10s 

♦ (b) 
^ J 

Fj6r(s) 

V. 

(A5) 

(AS) 

(A7) 


Transfer  Functions 


for 


By  solving  the  three  simultaneous  equations  (A5),  (A6),  and  (A7) 
0(s)  there  is  obtained 


^C0's5  + C-l's1*-  + C2's3  + C3's2  + C^'sj  0(s)  = 

+ Cg's^  + Cy's2  + Cg'sjbr(s)  (A8) 


where  the  Cn'  coefficients  are  defined  in  table  II.  Dividing  equa- 
tion (A8)  through  by  Cq'  and  using  the  unprlmed  coefficients  Cn  to 


designate 
and  C2  = 


the  ratio  of  the  primed  coefficient  to  Cq 
C2' 


/ c, 

( e • 8 • 7 C1  = ^T 


results  in  the  following  transfer  function  for  0 


0 

due  to 


a rudder-deflection  input  Bj.: 


p C^s^  + Cgs2  + CyS  + Cg 
6r  + c^s^  + CgS2  + C3S  + C 4 


(A9) 


Determined  in  a similar  manner,  the  transfer  function  for  0 is 

C9S2  + Cjqs  + 


£ = _ 

&r  + C}_s2  + C2s2  + C3S  + C4 


and  the  transfer  function  for  \|r  is 


JL  = C12s?  + C13s2  + Cl48  + C1‘5 

&r  sfs^  + C-j_s^  + C2s2  + CjS  + CjJ 


(A10) 


(All) 
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The  lateral-acceleration  transfer  function  may  be  developed  from  the 
sideslip  transfer  function  in  the  following  manner:  The  lateral  accelera- 

tion is  usually  defined  as 


ay  = v(Dp  + D#  - Kg0) 

Substituting  equation  (Al)  into  (A12)  yields 


(A12) 


ay  = VFlBr(t)  - VKlP 


(A13) 


Therefore, 


or 


P(s) 

&r(s ) 


VKX 


(Al4) 


fi  _ Cl6S^  + C17s^  + Cl8S2  + CjgS  + CgQ 
8r  s4  + c^s3  + C2S^  + CjS  + 


(A15) 


Modes 

The  characteristic  equation  of  lateral  motion  is  obtained  by 
expanding  the  determinant  of  the  coefficients  of  equations  (Al),  (A2), 
and  (A3)  and  is 

s^  + C-^s^  + c2s2  + CjS  + = 0 (Al6) 

This  equation  can  be  factored  as  follows: 

s(s  + Aj^s  + Ag^s  + I - 1tj)(s  + i + 1tj^=  0 (A17) 

where  the  root  s = 0 indicates  that  the  aircraft  is  insensitive  to 
azimuth,  the  root  A-j_  is  the  damping-in-roll  mode,  the  root  Ag  is 

the  spiral  mode,  and  the  root  -{■  + iq  is  the  Dutch  roll  mode,  an 
oscillatory  motion  composed  of  roll,  yaw,  and  sideslip.  The  modes  of 
this  airplane  are  indicated  in  table  X. 


t 


' i 


\ 
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An  approximate  expression  for  which  gives  excellent  results  is 


M = £itl 

2 n_  1 


(A18) 


Frequency  Response 

The  frequency  response  which  is  the  steady-state  response  to  a 
sinusoidal  input  consists  of  an  amplitude -ratio  and  phase-angle  relation 
between  input  and  response  for  various  values  of  co.  It  is  determined 
by  substituting  s = ico  into  equations  (A9)>  (A10),  (All),  and  (AI5), 
respectively,  for  3,  0,  \|r,  and  ay  and  is  a complex  expression.  The 

amplitude  ratio  is  the  square  root  of  the  sum  of  the  squares  of  the  real 

and  imaginary  parts  of  the  complex  expression  and  the  phase  angle  4 is 

the  arc  tangent  of  the  ratio  of  the  imaginary  part  to  the  real  part. 

The  amplitude  ratio  for  sideslip  angle  3 is 


p _ (c8  ~ c6a^)2  +(Cj“>  ~ Cjgyj2 

§r  ^ (to*  - C2a£  + C4)  2 + jcjcu  - C1Co5) 

and  the  phase-angle  relation  between  3 and  6_  is 


(A19) 


fCjw  - C^aA 

= tan--1-  — — - tan 

**  \°8  - 

The  amplitude  ratio  for  roll  0 is 


1 / Cjco  - Cjyy> 
lo)^  - Cgu£  + Cj^ 


(A20) 


jjf  = (cll  - Cep?)2  +(CiQa>)2 

&r  \ - Cga£  + C^2  + (cyn  - C ±<J> 

and  the  phase -angle  relation  between  j l and  6r  is 


(A21) 


4 = tan-1  ( — — — — \ - tan" 

^ lCll  “ C9°7 


Cyu  - C^uP 


o>k  - C2o^  + Ck 


(A22) 
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The  amplitude  ratio  for  yaw  is 


Jl 



5r 

'\ 

^15  " ^13°^)  ^ + - C-|  ?txr*j  ^ 


and  the  phase -angle  relation  between  ^ and  5r  is 


4 1 = tan-1  (CjAy --£12^1 . tan-1  r.  - ^ + M 


'15  “ u13a 


C^CD  - Cj(j£ 


The  amplitude  ratio  for  lateral  acceleration  ay  is 


h 

§r 

= \ 

|Ci6afr  - Cl8^  + C2Qj 2 + {C19a>  - C17aJj 


and  the  phase -angle  relation  between  and  6r  is 


4a  = tan 


\ f 

\ *»  t 

- - tan" 

16^  - Cl8a^  + C20/  r " C2U^  + cv 


-1 1 C19a>  - C17c0?  \ +an-i  / c3a)  - Cj^ 


(A23) 


(A24) 


(A25) 


(A26) 
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APPENDIX  B 

DETERMINATION  OF  LATERAL-STABILITY  DERIVATIVES  BY  USE  OF 
FREQUENCY  RESPONSES  DERIVED  FROM  CALCULATED  TRANSIENT 
RESPONSES  TO  AILERON  DEFLECTION 

Transient  responses  in  roll,  yaw,  and  sideslip  to  a square-pulse 
aileron  deflection  were  calculated  for  the  airplane  whose  mass  and  aero- 
dynamic characteristics  are  given  in  table  XI.  These  transients  and 
assumed  aileron  time  histories  are  presented  in  figure  J. 

The  frequency  responses  D^/&a,  Ety/sa,  and  3 /&a  were  obtained 

from  analysis  of  the  calculated  transient  responses  D <f>,  D|r,  and  3 to 
the  assumed  aileron  deflection  by  use  of  the  equation 


1 — St  Ti 


(t)dt 


e",stq(t)dt 


G(s) 


The  Laplace  transforms  of  the  responses  were  obtained  by  fitting  various 
functions,  such  as  polynomials  or  trigonometric  functions,  to  finite 
sections  of  the  respective  output  curves  and  thus  performing  analytically 
the  required  integrations  from  0 to  » in  finite  time  intervals  as 
described  in  reference  8.  These  transforms  were  then  evaluated  for  a 
number  of  values  of  io>  (where  s = icu)  and  frequency-response  data 
were  obtained.  These  data  have  been  plotted  in  figure  8 where  they  are 
represented  by  the  triangular  test  points. 

For  purposes  of  comparison,  the  frequency  responses  for  the  airplane 
were  calculated  from  the  basic  mass  and  aerodynamic  characteristics  and 
are  also  presented  in  figure  8 in  the  form  of  curves.  The  agreement 
between  these  two  sets  of  data  is  seen  to  be  excellent. 

In  addition,  the  frequency  responses  were  obtained  from  the  calcu- 
lated transients  by  use  of  HM  equipment.  The  Fourier  transform  of  each 
transient  was  obtained  from  a numerical  evaluation,  by  means  of  the 
IBM  equipment,  of  the  integral 


-ioJtp(t)dt 


' i 


l 

\ 


- . 


i 


* 
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between  the  limits  of  t = 0 and  t = 6.9  seconds  to  which  was  added 
an  end  correction  based  on  an  analytical  expression  of  the  transients 
from  t = 6.9  seconds  to  infinity.  A time  interval  At  = 0.05  second 
was  used  in  the  numerical  integrations.  The  frequency  responses  thus 
obtained  are  shown  in  figure  8 as  square  points  and  appear  to  be  in 
excellent  agreement  with  the  frequency  responses  calculated  from  the 
basic  mass  and  aerodynamic  characteristics  and  with  those  derived  from 
the  transients  by  the  curve -fitting  method. 


By  use  of  equations  (27),  (31)>  (33)>  (35)>  and  (36)  of  the  body  of 
the  paper  and  the  frequency-response  data  derived  from  the  transients  by 
both  methods,  the  coefficients  Ky,  Kj,  K^,  Kg,  G2,  Ky,  K^,  and  Ky0 

were  obtained,  and  hence  the  stability  derivatives  Cy^,  C^,  C^,  C jy, 

CZ&  , Cnp,  C^,  and  C^.  These  results  are  presented  in  table  XII 

along  with  the  values  assumed  for  the  calculation  of  the  transient  motions. 
Results  are  also  presented  for  the  case  where  equations  (l),  (2),  and  (3) 
of  the  body  of  the  paper  were  not  separated  into  real  and  imaginary  parts 
but  were  treated  as  vector  equations  and  the  method  of  least  squares 
applied  in  a manner  similar  to  that  described  in  reference  4.  The  amount 
of  work  required  in  the  least-squaring  process  is  less  when  these  equa- 
tions are  treated  as  vector  equations,  but  a limitation  is  introduced  in 
that  a parameter  such  as  Kg  which  is  difficult  to  extract  with  good 

accuracy  cannot  be  eliminated  as  it  was  in  the  body  of  the  paper  when  the 
equations  were  separated. 


A comparison  of  the  values  of  the  derivatives  based  on  frequency 
responses  obtained  by  the  curve-fitting  method  and  obtained  from  the  two 
methods  of  treating  the  equations  of  motion  indicates  very  good  agreement 
for  all  the  parameters  except  C7  . The  poor  agreement  for  C,  is  felt 

not  to  be  significant  since  C generally  makes  a negligible  contribu- 
tion to  the  lateral  motions  of  an  airplane  and  is  important  only  in  the 
lightly  damped  spiral  mode  of  motion  which  is  very  difficult  to  define 
accurately  from  a transient  response. 


A similar  comparison  of  the  two  sets  of  values  based  on  frequency 
responses  obtained  with  the  IBM  equipment  shows  good  agreement  except  for 


the  derivatives  C7  , C_  , 

4r  Uy 

initially  assumed  values  of 


The  poor  agreement  between  the 
, and  Cv  and  those  calculated 


-J  y V_/_  y CULiVU  Uy  auu  WUUOW  VCkLU  uxavcu 

iT  V 

from  the  frequency  responses  obtained  with  the  IBM  equipment  is  probably 
due  partly  to  the  fact  that  data  were  not  obtained  for  frequencies  near 
the  natural  frequency  of  the  airplane  with  the  particular  routine  used 
for  the  IBM  equipment  and  partly  to  small  inaccuracies  that  exist  in  the 
data  at  the  frequencies  that  were  included.  These  small  inaccuracies 


V /■> 
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undoubtedly  could  have  been  reduced  somewhat  by  the  use  of  a smaller  time 
interval  in  the  numerical  integrations  performed  by  the  IBM  equipment. 

The  use  of  the  frequency-response  data  which  were  obtained  from  the 
transient  responses  by  the  curve-fitting  method  and  which  included  fre- 
quencies near  the  natural  frequency  of  the  airplane  resulted  in  good 
agreement  for  these  derivatives.  Therefore,  it  appears  that  further 
investigation  into  the  effects  of  the  choice  of  the  frequency  range  to 
be  used  in  extracting  stability  derivatives  from  frequency-response  data 
is  warranted. 
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TABLE  I.-  DEFINITIONS  OF  THE  COEFFICIENTS  OF  THE 


EQUATIONS  OF  LATERAL  MOTION  OF  AIRPLANE 


Stability -derivative  coefficients : Control -effectiveness  coefficients : 


% =%(-£) 

k2  * * e/v 


K3  = % 


r. a — 

L 2(kx/b)2  T2_ 


°Zp  _ 

% * ixz/xx 


Fi  - V b 


f2  = c 


16'2m2t2 


f3  - 


r 2 M 2 t2 


Gp  — G n ■■■  * 

* 2 (kx/b)  T 


“6  = °l, 


4t  (kx/b) 


Kq  = Ixz/lZ 


“9  -%S7 


K10  = cnr 


'(kz/b) 
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TABLE  II.-  COEFFICIENTS  OF  TRANSFER  FUNCTIONS 


The  equations  for  the  transfer  functions  are  given  in  appendix  A; 
the  coefficients  Cn  are  obtained  from  the  coefficients  Cn' 

by  dividing  by  C0'  |e.g.,  Cx  = 

C0*  - 1 - 8 

Ci'  =K1  + Kk  + K10  - KjKj  - K^Kq  - K^Kg 

C2'  = K?  + KXK4  + KXK10  + K4K10  - KjKg  - - KjKjKj  - K^Kg 

C3'  = KgK^  + - K5Kg  + K1K4K10  - KjKgl^  - 

<V  - K2K3k10  - W<7 

C5 ' = Fl(x  - k5kq) 

c6'  = FjK4  + FlKl0  - K8F2  - F3  - F 1KjS9  - FjKgKg 
C?  ’ = FxK4K10  - F2K$  + FjKgK^  - FjK4  + FgKg  - FjKgl^ 

C8'  « + f2k10k2 

c9'  = f2  + F^ 

Cio'  = FgKi  + F jK^Krj  - FxK3  + FjKg  + FjKjK^  + F^q 
Cil’  - F2KjK10  + FjKgCy  + F3K3  + F^  - F^Kio  + Fj^Kg 
Ci2‘  = F3  + KgF2 

C1  y = F3K4  + F5Kx  - F1K5Kg  + F^  + F^  + F^Kg 
Cl4’  = F3KxK4  - FjKjKg  + F1K4K7  + 

^15 ' = FgKglC^  + FjKgKj 

Cl6’  = Co'VFx 

C17'  = C^VFj.  - VKjCy 

Ci8'  = VFjCg’  - VKjCg’ 

Ci9'  = VF^'  - VKj<y 
C20'  = VFxC4'  - VKjCq’ 


' i 
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TABLE  III.-  AIRPLANE  PARAMETERS  USED  IN  EXAMPLE  GIVEN  IN  BODY  OF  PAPER 

(a)  Known 

M = 0.8  (evaluated  for  an  altitude  of  10,000  ft) 
b = 22.6  ft 
S = 130.0  sq  ft 
p = O.OOI756  slug/cu  ft 
V=  861.74  ft/sec 
Ix  = 2062  slug-ft2 

Iz  = 13,298  slug -ft2 
m = ^ = 295.03  slugs 

Ixz  = 157  slug-ft2 

cnp  = 0 

(b)  Computed 

T " Jv  “ x‘5  sec 
- ”-2 


K = = iffi  = 0.07614 

Ix 

K8  = ^ = 0.011806 

XZ 

K2  = | = 0.037^  per  sec 
- ■ • 


1 
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TABLE  IV.-  FREQUENCY  -RESPONSE  DATA  FOR  EXAMPLE 
GIVEN  IN  BODY  OF  PAPER 

(a)  Sideslip,  3 


O), 

radians 

sec 

L 

Sr 

®fV 

deg 

A3  = 

&r 

cos  $>3&r 

Bp  = 

P 

&r 

sin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0.536 
• 574 
.650 

.797 

1.116 

2.105 

4.338 

1.368 

.715 

.464 

-0.1 

-1.1 

-2.3 

-4.1 

-8.2 

-20.8 

-109.7 

-158.1 

-166.1 

-169.1 

0.536179 

.574106 

.649929 

.794967 

1.104368 

1.967417 

-1.462174 

-1.268821 

-.693781 

-.455533 

-0.001562 

-.011025 

-.026083 

-.056988 

-.159106 

-.747357 

-4.083763 

-.510100 

-.171676 

-.087720 

(b)  Roll, 


CD, 

radians 

sec 

i 

5r 

V 

deg 

h 

l 

sr 

cos  ^ 

B0  = 

1 

Sr 

sin  4^ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

8.761 

4.624 

3.432 

3.111 

3.459 

5.446 

9.719 

2.729 

1.298 

.780 

79-0 

67.7 

57-4 

47.6 

37-1 

I8.9 

-74.6 

-126.9 

-138.2 

-143.9 

1.676642 

1.754686 

1.850156 

2.096761 

2.760823 

5.153317 

2.578202 

-1.638908 

-.966715 

-.630341 

8.598903 

4.278325 

2.891147 

2.297575 

2.084203 

1.762687 

-9.370572 

-2.181509 

-.865863 

-.459096 

i 


1 

1 


% 
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TABLE  IV.-  FREQUENCY -RESPONSE  BATA  FOR  EXAMPLE 
GIVEN  IN  BODY  OF  PAPER  - Concluded 


(c)  Yaw,  \|f 


(d)  Lateral  acceleration,  ay 
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TABL*  V.-  IBDiraSD  FORM  OF  EQUATIONS  USED 
M EXAMPLE  OIVDi  IN  SOW  OF  PAPER 

- *.  * 

“El  ' (K]*s  - KgA^  - uiBp  - a&^  ^ 


(27»> 

(aa.) 

(29») 


(50*) 


S^iHJ  ^NjW, 

V 

* «« 

(A  - ^)j(.Be)j 

KiO 

rJ,3  ‘ ^ ('Vs  * *6“%  * *10“**  " ^\)  3 


j^Ev.  ***)]• 

^(“V.  - “V»)j(*eN  * V*)j 

^(^(WVs), 

lt5 

£ 

j-h 

(H** + *»*»)  j (“*A  - «*#S)  j 

£((-A--*)J* 

*fc 

_ 

£kt 

Fp 

* Vt)j 

j pv.  * °V*  - ^|tf)j  HA  ' “V*)  j } 

^ (f2  * ^ - K^A,  ♦ *„<*,  - KjAp)  J f-V)  ^ 

^(v»  ♦ V*d 


(»*•) 


(55») 


(56.) 
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TABLE  VI.-  LEAST-SQUARES  SOLUTION  OF  EQUATIONS  FOR  DATA 
OF  EXAMPLE  GIVEN  IN  BODY  OF  PAPER 


(a)  K]_,  equation  (27a) 


£(- 


aAp  - aAy  + KgB^Bp  - 7*378 


^(B?)1 2  . 17.562 


17.562K].  = 7.378 
Kx  = 0.420 


(b)  F-^,  equation  (28a) 

(KlAp  - KgA^  - oBp  - oB^  = 1.06l 


10Fl  = 1.061 

F1  = 0.106l 


1 

* 


V 
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13,042,000%  = 5,568,000 
% = 0.427 

(d)  % from  lateral -acceleration  data,  equation  (30a) 

ZZ(K1^  + ^)=  1.035059 

OJ=l  v v ' 

10F1  = 1.035 

Fx  = 0.1035 

\ 
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TABLE  VI.-  LEAST-SQUARES  SOLUTION  OF  EQUATIONS  FOR  DATA 
OF  EXAMPLE  GIVEN  IN  BOEV  OF  PAPER  - Continued 


(e)  Ky  and  Kj_q,  equation  (31a) 

AfL  , x 2 

(-  bb)  * = 17.562397 

oj=1  ' ' 

B^aA^  = 44.47131 


473.928987 


Bp)  - 856.325838 
(o)%  - K8u)2B^(oA|i)  = 2356.791416 


Kq  = 0.011806  Ky  = 0 
17.561c,  + 44.47K10  = 856.3 
44.471C,  + 473-95K10  = 2356.8 
Kj  = 47.443 
K1q  = 0.5217 


(f)  Fj,  equation  (32a) 

IOF3  = KyAp  + Kgo^  - - a^A^ ) = -251. 87 


Fj  = -25.187 


' i 

l 

: ' \ 
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TABLE  VI.-  LEAST-SQUARES  SOLUTION  OF  EQUATIONS  FOR  DATA 
OF  EXAMPLE  GIVEN  IN  BODY  OF  PAPER  - Continued 


(g)  K4,  and  F2,  equation  (34a) 

+ A^U2  = 371.762380 

a^l  ' 


10 

2Z(BpB^  + ApA^jo^B^  - oB^)=  -3617.677695 


10 


o>=l 

10 


(BpB^  + ApA^(-  A^=  23.184790 


- oB^)2  = 35234.438878 


Air ) = -209.219323 


^ (-  aJ)  2 = 11.122853 

o>=l  V 7 

+ “V  - %»2|||^  (3  By  + ApA^)  = 33187.73450 9 

" uBA)=  -522344.219155 

^B0By  + a&ySr  " ^)("  \)  = 2422.633922 


371.76238K5  - 3617.67769K4  + 23 • l8479F2  = 35187.75451 
-3617.67769KJ  + 35234.43888^  - 209.21932F2  = -322344.21915 
23.18479K5  - 209.219321^  + 11.12285F2  = 2422.633922 
Kj  = 138.272  K4  = 5.212  Fg  = 27.636 


( 


% 
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TABLE  VI.-  LEAST -SQUARES  SOLUTION  OF  EQUATIONS  FOR  DATA 
OF  EXAMPLE  GIVEN  IN  BODY  OF  PAPER  - Concluded 

(h)  Kg,  equation  (35a) 

Z.M2*  811.55627 

cupI  ' 

^(a^)(F2  + o)2A^  - K^2^  + Kka£p  - K?Ap)  = 235-117592 

811. 6Kg  = 235.1 

Kg  = 0.290 

(i)  Alternate  means  of  obtaining  Kg,  equation  (36a) 

trK)2  = 473.928983 

+ K4oA^  + K^B^  - 137-77^056 

473. 9Kg  = 137.77 
Kg  = 0.291 
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TABLE  VII.-  COMPARISON  OF  COMPUTED  AND  KNOWN  VALUES  OF  COEFFICIENTS  OF 
EQUATIONS  OF  LATERAL  MOTION  FOR  EXAMPLE  GIVEN  IN  BODY  OF  PAPER 


Coefficient 

Computed  from  data 

Known  value 

K1 

J ®0. 42054 
l b0.4269 

| 0.427 

*3 

138.272145 

138.245 

*4 

5.212568 

5.21 

*6 

O.290 

0.3017 

k7 

47.438625 

47.41 

*10 

0.521457 

O.5272 

*1 

J *0.10609 
[ b0.1035 

| 0.104 

*2 

27.636187 

27.65 

f3 

-25.187091 

-25.22 

aFrom  sideslip  data. 

bFrom  lateral -acceleration  data. 


TABLE  VIII.-  COMPARISON  OF  COMPUTED  AND  KNOWN  STABILITY 
DERIVATIVES  FOR  EXAMPLE  GIVEN  IN  BODY  OF  PAPER 


Stability 

derivative 

Computed 

Known 

% 

-1.259 

-1.28 

C*P 

-0.1490 

-0.149 

C*P 

-0.4280 

-0.428 

Czr 

0.02389 

0.0248 

cn$ 

0.3292 

0.329 

Cnr 

-0.2732 

-0.279 

O.3I83 

0.312 

0.02979 

0.0298 

C% 

-0.1748 

-0.175 

\ 


i 


1 


% 


r 


■ & 


\ 


v 


! 
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TABLE  IX.-  COMPARISON  OF  COMPUTED  AND  KNOWN  COEFFICIENTS  OF 
TRANSFER  FUNCTIONS  FOR  EXAMPLE  GIVEN  IN  BODY  OF  PAPER 


Coefficient 

Computed  from  data 
(least-squares  value) 

Known 

Co' 

0.999101 

O.999IOI 

Cl’ 

6.150763 

6.160254 

C2* 

50.934260 

50.972855 

C3' 

253.^3361 

253.2143 

C4‘ 

2.193801 

2.190861 

c5* 

0.105995 

O.IO3907 

C6’ 

25.468778 

25.48986 

C 7' 

132.544801 

132.64415 

C8’ 

O.267223 

0.260611 

<V 

25.718442 

25.729749 

cio' 

3.636361 

3.952748 

Cn' 

-2174.851610 

-2178.769199 

Cl2* 

-24.860818 

-24.893564 

C13' 

-136.884812 

-137.264852 

ClV 

-28.978823 

-30.417543 

C15' 

-81.654072 

-81.36940 

cl6' 

89.IH765 

89.62096 

C17- 

509.187633 

514.350717 

cl8' 

-4822.026985 

-4815.419162 

C19* 

-26127.620554 

-26138.324577 

c2o' 

97.287337 

100.542037 

TABLE  X.-  COMPUTED  LATERAL  MODES  FOR 

EXAMPLE  GIVEN  IN  BODY  OF  PAPER 

[Characteristic  equation  is  equation  (Al6)  of  appendix  a] 

Damping-in-roll  mode,  -5*393 

Spiral  mode,  Xg * -0. 008668 

Dutch  roll  mode,  -£  ± irj -O.381  + i(6.84) 


J 


i 


\ 

\ 


\ 


/ 

1 


V 
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TABLE  XI.-  PARAMETERS  OF  AIRPLANE  USED  AS 
EXAMPLE  IN  APPENDIX  B 

Altitude,  ft  30,000 

Wing  loading,  lb/sq  ft 65 

V,  ft/sec 797 

b,  ft 28 

CL 0.23 

Mb 80. .7 

Kx2  0.0097 

KZ2 0.0513 

Kxz 0.00145 

per  radian -0.40 

C2r  per  radian 0.08 

Cn^  per  radian -0.02 

Cnr  per  radian -0.40 

Cv.  per  radian -1.0 

p 

Cn„  per  radian 0.25 

p 

C7  per  radian -0.126 

C per  radian -0.11 
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Figure  3«-  Frequency  response  in  sideslip  due  to  rudder-deflection  input 
for  rigid,  high-speed  airplane  of  example  given  in  "body  of  paper. 


Figure  4.-  Frequency  response  in  roll  due  to  rudder-deflection  input 
for  rigid,  high-speed  airplane  of  example  given  in  body  of  paper. 


Figure  5.-  Frequency  response  in  yaw  due  to  rudder-deflection  input  for 
rigid,  high-speed  airplane  of  example  given  in  body  of  paper. 


Figure  7.-  Calculation  of  transient  responses  to  a square-pulse  aileron 
deflection  for  the  airplane  defined  in  table  XI. 


Figure  7 — Continued. 
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Figure  8.-  Continued. 


Calculated  from  mass  and  aerodynamic 


Figure  8.-  Continued 
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